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Abstract 

We study the finite temperature electroweak transition with non-perturbative lattice 
Monte Carlo simulations. We find that it is of first order, at least for Higgs masses 
up to 80 GeV. The critical temperature of the phase transition is found to be smaller 
than that determined by a 2-loop renormalization group improved effective potential. 
The jump of the order parameter at the critical temperature is considerably larger 
than the perturbative value. By comparing lattice data and perturbation theory, we 
demonstrate that the latter, for the computation of the vacuum expectation value of 
the Higgs field v(T) in the broken phase at given temperature, converges quite well, 
provided v(T)/T > 1. An upper bound on the Higgs mass necessary for electroweak 
baryogenesis in the light of the lattice data is briefly discussed. 



CERN-TH.7244/94 
May 1994 

Martially supported by a CEC Science program (SCI-CT91-0729). 

2 On leave of absence from the Institute for Nuclear Research of the Russian Academy of Sciences, 
Moscow 117312, Russia. 



1 Introduction 



The high-temperature limit of the 4-dimensional electroweak theory near the phase 
transition and for sufficiently large Higgs masses can be described by an effective 3- 
dimensional gauge-Higgs theory (for a discussion of high-temperature phase transitions 
see IJ, for a general consideration of dimensional reduction see [§], for 3d EW theory 
see P, f|, U). The effective theory is strongly coupled in the symmetric phase || |7j; 
in particular, pure SU(2) gauge theory in 3d is confining. This fact makes a reliable 
entirely perturbative computation of almost all characteristics of the phase transition 
impossible. For example, the computation of the critical temperature T c requires the 
comparison of the vacuum energy of the symmetric and broken phases, and the former 
cannot be estimated on perturbative grounds. The bubble nucleation temperature T* 
in cosmology is related to the structure of the effective action for small scalar fields, 
which is precisely the place where infrared divergences are most severe. The expectation 
value of the Higgs field v(T*) at the temperature T*, necessary for the estimate of the 
fermion number non-conservation rate f|, cannot be found perturbatively, just because 
T* is non-perturbative. At the same time, the determination of the vacuum expectation 
value of the Higgs field in the broken phase at a given temperature T by perturbation 
theory does not require any information on the symmetric phase and may be good 
enough provided that the ratio v(T)/T is sufficiently large. Therefore, the study of 
the nature of the electroweak phase transition requires non-perturbative methods. The 
most straightforward ones are lattice Monte-Carlo simulations. 

First lattice Monte Carlo results on the electroweak phase transition using a 3d 
effective theory have already been given in (for 4d simulations see ||). We improve 
here the analysis of 0] in two ways: by going from the 1-loop to the 2-loop level in 
the discussion of the 3d effective theory and its effective potential, and by extending 
considerably the numerical calculations. We find that the EW phase transition is of 
the first order at m# ~ mw- The critical temperature we determined is smaller than 
that following from the perturbative analysis, while the jump of the order parameter is 
larger. At the same time, 2-loop effective potential gives an adequate description of the 
order parameter in the broken phase for v(T)/T > 1. These results are in a qualitative 
agreement with a picture of the electroweak phase transition suggested in [|nj . Given 
the information provided by lattice simulations we comment shortly on the question of 
the Higgs mass necessary for electroweak baryogenesis |TT], [TTJ. 

The advantages of considering the effective 3d vacuum gauge-Higgs system in com- 
parison with the full 4d theory at non-zero temperatures have already been discussed 
in M (from the lattice point of view) and in || (from the perturbation theory point of 
view). In the present paper, for reasons of simplicity, the U(l) factor of the complete 
electroweak theory is disregarded. The fermionic contributions are omitted as wellQ. 

The paper is organized as follows. In Section 2 we briefly summarize the results 

3 The most important fermionic contribution is that of the top quark to the effective scalar mass 
at high temperatures. It does change the absolute value of the critical temperature, decreasing it, 
but changes only marginally dimcnsionless ratios such as v(T c )/T c , most important for cosmological 
applications. 
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from the perturbative analysis. The lattice action is given in Section 3, as well as a 
discussion on relating lattice numbers to physical quantities. Numerical results are 
presented in the Section 4. Section 5 is discussion. 



2 Continuum Lagrangian 

The continuum Lagrangian of the effective 3d theory under consideration is 
S cS = J dPx^FZFS + i(AA ) a (AA ) a + (A0) f (A0) + 

\mlAlA% + \\a{A%A%? + mfrV + H<I>W + MS^V}- (1) 

It contains gauge fields together with the Higgs doublet and a scalar triplet (the former 
time component of the gauge field). Here all bosonic fields have the canonical dimension 
[GeV] 2 and 3d gauge and scalar couplings g%, A 3 , A^ and h 3 have dimension [GeV]. 
The relation between 4d and 3d coupling constants and masses on the 1- and 2-loop 
levels has been discussed in ||. For all numerical and analytical work we use the 
following simplified relations: 

<7 = ~, 9l = 9 2 T, \ 3 = lg 2 3 ^, h = ±gl X A = 0. (2) 

Due to the super-renormalizability of the 3d theory (only a finite number of diagrams 
are divergent) these relations are scale-independent^. The Debye screening mass and 
the effective Higgs mass are || 

ml = \g 2 T\ (3) 
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The numerical constant c = —0.348725, the parameter 113 is the scale of the MS scheme. 
The Debye screeening mass mo is a 3d renormalization group invariant, but the effective 
Higgs mass runs with the normalization point due to logarithmic divergences on the 
2-loop level. The 2-loop coefficient /2m is given by 

f2m = ^3 + 9A 3 S 3 2 " 12Ai (5) 

The 2-loop renormalization group improved effective potential for this theory has 
been computed in ||. (The computation of the high-temperature limit of the 2-loop 
effective potential in 4d has been performed in fT2| ). We do not write the expression 
here for lack of space and just mention that, for m H = 80 GeV and m w = 80.6 GeV, 
the critical temperature and the expectation value of 4> at T c in the broken phase 
computed from it are T c = 173.3 GeV and v(T c ) = 81 GeV (v(T c )/T c = 0.47). 

4 The actual value of the coupling is rather small, 17 g 4 T / (48tt 2 ) . We included this coupling in 
the lattice simulations just to make sure that the lattice action is bounded from below for constant 
fields, and checked that the results are independent of Xa provided that it is small enough. 
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3 The lattice action 

Going over to a 2x2 matrix representation $ = (0o + io^i) / y/2 of the doublet scalar 
field, discretising and scaling the continuum fields by 



igaA c 



Ac 



a 2 ' 



(6) 



the lattice action corresponding to the continuum action in eq. ([T]) becomes, in standard 
notation, 

x i<j z 

+k^^[Tr J 4„(x)Kr 1 (x)io(x + l )[/ 1 (x) -Tr^(x)] + 
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(7) 



where E = 0.252731. 

All the three lattice coupling constants are given in terms of g, T and win by the 
following equations, which directly follow from the discretisation procedure: 
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(9) 



The tree relation between the mass parameter m 3 and the lattice variables is given by 



rrin 



2(l-20 R -3/3 H ) 
(3 H a? 



(10) 



The 1-loop counterterm (~ l/o) removing divergences from the Higgs mass has been 
determined in [ffl, and the 2-loop counterterm ~ log(a) has been estimated in ||. For 
present lattices the linear counterterm, included in [|J, is clearly the dominant one, 
but the accuracy of the calculations is already such that the effects of the logarithmic 
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counterterm are clearly seen. No other counterterms exist since the 3d theory is super- 
r enormalizable . 

The knowledge of these counterterms allows one to relate lattice parameters to the 
temperature: 
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Because of the super-renormalizability of the 3d theory, this relation is exact in the 
continuum limit. On the 1-loop level the analogous relation has been derived in M. 
The determination of the (m^/m^ dependent) constant rj requires a computation of 
the 2-loop effective potential of the theory by the lattice perturbation theory and a 
comparison of the result with the corresponding expression in the MS scheme. Due to 
the complexity of this computation, it is not attempted here. Instead, we determine 
the constant rj by Monte-Carlo methods. 



4 Lattice results 

The continuum limit corresponds to Pg — * oo, N/Pg oo, with iV being the lattice 
size. Too large a (3q makes the system hard to simulate and we have used fta = 
12, 20, 32. Our computer resources have permitted us to perform runs on lattices of 
sizes 8 3 ... 32 3 . The confining property of the symmetric phase must not be lost, and 
this requires 

< 1.468JV, (12) 



at least for pure 3d SU(2) gauge theory |13[. The use of different values of Pg and N 



is important to study the scaling and consistency of the results. 

Simulations have also been carried out for the Higgs masses mj? = 35, 60, 70, 80 
and 90 GeV. We will discuss here in some detail only those simulations with 80 GeV 
Higgs (other Higgs masses are discussed in Qj). 

The first question is the very existence of the first-order phase transition. To an- 
swer it, we searched for a two-peak signal in the distributions of the different order 
parameters (such as length of the Higgs field R 2 , the average of the link operator 
L = TrV\x)Ui{x)V{x + 1), where $ = RV, etc.) A typical picture at p H = 0.347710, 
Pg = 12 and V = 24 3 is shown in Fig. 1. The two-peak structure is clearly seen. At 
the same time, due to the fact that the lattice volume is finite, the width of any of the 
peaks is comparable with the distance between them. For comparison, we present in 
Fig. 2 the distribution of R 2 for a Higgs mass ran = 35 GeV, where finite-size effects 
are less important. 

The determination of the critical value of P c u in the continuum limit requires the 
analysis of finite-size effects. In principle, a number of different methods can be used 
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for its determination: 

(i) Equal-area signal. The value of (3 C H is determined as the value of (3h for which the 
areas [15[] under the two peaks in the distribution of some order parameter (actually, 
it does not matter which order parameter is chosen) are equal. 

(ii) One looks for a maximum of the heat capacity (action susceptibility), considered 

m 



as a function of Ph 
(in) The same as in 



X5|. 



, but for L-susceptibility Jl6|, |TT|| . 
All three methods must give the same results in the continuum limit. For finite 
lattices and for finite Pa the results are different, and the continuum limit should be 
found by extrapolation. 

The first method works well for small Higgs masses, such as 35 GeV, since the two 
peaks are well separated, but it cannot be applied with good accuracy to the study 
of the phase transition with heavy enough Higgses, at least for the lattice sizes we 
used. So, we used (ii) and (iii) to determine P° H . To find the maximum values, we 
used the Ferrenberg-Swendsen multihistogram method |T7]], which yields Monte Carlo 
observables as continuous functions in /3jj around the actual Ph -values used in the 
simulations. The error analysis was performed with the jackknife method. 

The results of the determination of the critical P H are shown in Fig. 3 (Pa = 12) and 
Fig. 4 (p G = 20) as a function of 1/iV 3 , with N = 8, 12, 16, 20, 24 and 32 (without 
the point N = 20 for j3 G = 20). One can see that both methods give consistent results 
for sufficiently large volumes of the lattice. 

The extrapolation of these results to infinite volume requires some care. For /3a = 12 
eq. ([i~2|) gives N > 8. The (5a — 12 data for N > 12 can be fitted very well by a linear 
function (3 C H (V) = P c H (oo) + c/V, giving p c H (oo) = 0.347703(10). The power low fit 
P H (V) = P H (oc) + cV n for all lattice volumes gives a consistent result: (3 c H (oo) = 
0.347698(12) with the power n = 0.9 ± 0.1. 

For /3a — 20 finite-size effects are more important. From eq. ( |i~2"|) one gets for 
Pa = 20 N > 14, so that one to two points at the right of Fig. 4 cannot be trusted. 
Then, the linear fit with the use of the four largest volumes gives P H (oo) = 0.341721(9); 
with the three largest volumes, P H (oo) = 0.341710(ll);with the three largest volumes, 
P c H (oo) = 0.341668(21). 

Our final aim is the determination of the critical temperature. To fix it, we must 
know the parameter rj appearing in eq. (P"T|). Since there are no analytic computations 
of this parameter with the help of lattice perturbation theory (even if they are possible), 
we determine it by comparison of the results of lattice simulations with perturbation 
theory deep in the broken phase, where it should work well. We choose the order 
parameter R 2 for this purpose. In order to relate it to continuum physics one has to 
subtract from it linear (~ 1/a) and logarithmic (~ log(a)) divergent terms. It can be 
shown that (R 2 ) is related to the continuum effective potential through 
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(13) 



where the expression in square brackets contains terms that are divergent in the contin- 
uum limit, fi is an arbitrary normalization point (the total expression is /x-independent 
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due to the \i dependence of the effective potential V e ff), fj is unknown constant. The 
derivative of the effective potential is taken at in the minimum of the broken phase. 
In the leading approximation, ^^f\v(T) = \v{T) 2 . The two-loop counter-term con- 
tribution to expression ( |I3| ) is numerically suppressed, so that uncertainties in the 
determination of rj coming from an absence of information about fj are small provided 
fj is not too large (say, fj < 2). 

The evolution of R 2 distributions in the broken phase with Ph is shown in Fig. 5. 
In Fig. 6 we show (R 2 ) as a function of Ph both as given by lattice data and by 2-loop 
perturbative computation with parameter rj = 0.54 (Pg = 12) and 77 = 0.12 (Pg = 20). 
The Pg — 12 lattice data corresponds to the variation of v(T)/T from 1.3 to 2, and 
Pg — 20 data to v(T)/T = 1.7-3.9. One can see that the use of 2-loop effective 
potential for the determination of the vacuum expectation value for the Higgs field in 
the broken phase at given temperature works quite well (numerically < 1-2%). 

This means that the dimensionless expansion parameter ~ g 2 T/mw(T) is sufficiently 
small at these temperatures. In other words, non-perturbative effects and higher order 
terms can be essential numerically only for v(T)/T < 1. 

At the same time, the systematic difference between the lattice data and 2-loop 
continuum predictions is clearly visible (it is, though, very small, less than 1% in 
(R 2 ))- In Fig. 7 we show the the plot of Rf att — R 2 heor as a function of Ph and with 
the same 77 choices as in Fig. 6 (the best fits). We suspect that the main source of 
the difference comes from finite size and finite spacing effects, which are difficult to 
compute analytically. As a rough estimate of parameter rj we take 77 = 0.3 ± 0.5, the 
central value being the average of Pg = 12 and Pg = 20 best fit values, while the error 
estimate is the difference between them. 

With the value of rj determined, the observed values of P C H can be converted to 
results for the critical temperature of the electroweak phase transition using eq. (pTT]) . 
From the Pg = 12 data we get T c = 162.1 ± 2.6 GeV and from the Pg = 20 data 
quite consistent values 160.3 ±2.6, 160.9 ±2.8 and 163.4 ±3.0 GeV, depending on the 
extrapolation to infinite volume (using four, three, two largest volumes, respectively). 
One observes that the T c determined by lattice methods is clearly smaller than the 
value T c = 173.3 GeV determined by perturbative methods. 

Since we do not have a clear peak separation at the critical temperature (implying 
that the system contains, at this temperature, a large fraction of interface configura- 
tions interpolating between the broken and symmetric phases) we cannot unambigu- 
ously extract from the lattice data at T c the quantity (R 2 ) relevant for the broken phase 
only. However, if we take the position of the peak to the right of Fig. 1 as an estimate 
of {R 2 ) in the broken phase, then we see from eq. ( |T3"D that v(T c )/T c = 0.73 ±0.04 (the 
error estimate is based on the width of the peak). This value may be compared with 
the results computed in optimized 2-loop perturbation theory: v(T c )/T c = 0.47 and 
v(T)/T = 0.81 for, say, T = 163.6 GeV (which is taken as a "true" critical tempera- 
ture). We conclude that the true jump of the order parameter (0.73) is clearly larger 
than that (0.47) coming from perturbation theory. At the same time, the perturbative 
prediction of v (T) at T = 163.6 GeV is rather close to the true value (compare 0.73 
and 0.81). The difference may be due to the fact that the unambiguous extraction 
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of v(T)/T from the distribution of R 2 is not possible and that finite size effects are 
not small. Another possibility is that at v(T)/T ~ 0.7 deviations (perturbative or 
non-perturbative) from 2-loop potential predictions large as 10%. 



5 Discussion 

The lattice results indicate that the phase transition is more strongly first order than 
could be expected from perturbation theory. At the same time, a good convergence of 
the perturbation theory in the broken phase allows us to relate these deviations with 
the properties of the symmetric phase. At present we cannot extract from available 
lattice data the parameters relevant for cosmological applications. In particular, the 
region of metastability of the symmetric phase, together with the bubble nucleation 
temperature T* cannot be determined. The direct lattice determination of T* is hardly 
possible at all, so that some indirect methods are necessary. One of them is related to 
the measurement of the interface tension at the critical temperature, which allows one 
to compute the bubble nucleation rate at least in the vicinity of the critical temperature. 
Another obvious option is to study details of the phase transition for the heavy Higgs 
boson, since we are approaching the non-perturbative region from the side of the broken 
phase, increasing ran- We plan to return to these questions in future work. 



As has been discussed in |I0| , non-perturbative effects are likely to decrease the value 
of the potential in the vicinity of the origin. It seems that our lattice results are in a 
perfect qualitative agreement with this picture. Indeed, consider the plot of the 1-loop 
and 2-loop renormalization group improved effective potentials at the "true" critical 
temperature T* rue ~ 163.6 shown in Fig. 8 (note that the deviation of the 1-loop 
result from the 2-loop one is negligibly small thanks to the optimization procedure, 
0). Clearly, these potentials are completely wrong at small fields 0, since we know 
that we are at the critical temperature. Therefore, the effective potential must have a 
contribution making the broken and symmetric phases degenerate. According to Fig. 
8 this contribution at = equals 

V peit (v(T c )) - WO) ~ 0.03# 3 6 w m(a w Tf. (14) 

A possible non-perturbative contribution is shown in this picture by a dashed line 
(in notations of ref. |]10|, we have taken Ap = 0.36 in the non-perturbative part of 



the potential). If the specific model for the description of non-perturbative effects 
considered in |L(j] is correct, then with this value of the non-perturbative energy shift 
electroweak baryogenesis may be possible up to a Higgs mass of about 100 GeV. The 
exact determination of this bound, however, requires a much better understanding of 
the symmetric phase of the electroweak theory at high temperatures. 
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Figure 3: The critical values of (3h for mu = 80 GeV and (3g = 12 for different lattice 
sizes. The full points correspond to the maximum of L-susceptibility, and empty points 
to the maximum of the heat capacity. 
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Figure 4: The same as Fig. 3 but for (3 G = 20 
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P G = 12, m H = 80 GeV 




Figure 5: Evolution of the distribution of R 2 in the broken phase with (3 H for (3q = 12, 
with m H = 80 GeV. 
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Figure 6: (R 2 ) determined on the lattice (points) and by perturbation theory (solid 
lines). The upper curve corresponds to rj = 0.12, (5q = 20 and the lower one to 
r] = 0.54, p G = 12. 
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Figure 7: The difference between {R 2 } determined on the lattice and (R 2 ) found with 
the use of perturbation theory for rj = 0.12, (3 G = 20 and rj = 0.54, (3 G = 12. 
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Figure 8: The renormalization group improved effective potential. Solid line - 2- loop 
order, dashed line-l-loop order, short-dashed line - possible non-perturbative contri- 
bution. The x-axis is the 4d scalar field in GeV, the y-axis is the dimensionless effective 
potential 12g^ 6 V e //(0). 
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